Decoherence of a Superposition of Macroscopic Current States in a SQUID 
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We show that fundamental conservation laws mandate parameter-free mechanisms of decoherence 
of quantum oscillations of the superconducting current between opposite directions in a SQUID — 
emission of phonons and photons at the oscillation frequency. The corresponding rates are com- 
puted and compared with experimental findings. The decohering effects of external mechanical and 
magnetic noise are investigated. 
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I. INTRODUCTION 

The possibility of macroscopic quantum tunneling of 
the magnetic flux in SQUIDs was suggested in a seminal 
paper of Caldeira and LeggettEl and was subsequently 
demonstrated in experimentcrtl. Following these devel- 
opments the possibility of coherent quantum oscillations 
between macroscopic flux states was intensively studied 
by theorists (see, e.g. Ref. ^] and references therein). 
Modern interest to this problem was generated by the 
hope to build a superconducting (SC) qubit. The goal 
of preparing quantum superposition of macroscopic flux 
states in a SQUID had remained elusive- however, un- 
til recent experiments of Friedman et alE and van der 
Wal et all]. In these experiments the tunneling split- 
ting between the states corresponding to symmetric and 
antisymmetric quantum superpositions of clockwise and 
counterclockwise currents in a SQUID loop have been 
observed. These experiments were follosajed by similar 
measurements of more elaborate designsLrtl 

The most important question about tha. above exper- 
iments is the one of the decoherence timet! Any degree 
of freedom that interacts with the coherently oscillat- 
ing variable can be the source of decoherence. One can 
divide all such sources into two groups, avoidable and 
unavoidable. Examples of avoidable sources are, e.g., 
nuclear spins and non-thermal noise. In this paper we 
study generic mechanisms of decphpjence which are con- 
trolled by the conservation lawsE2lO. In application to 
SQUIDs, such mechanisms are unavoidable as they orig- 
inate from the fundamental symmetries of free space. To 
illustrate our point, consider, e.g., experiment of Ref. ^ 
in which the current, J^3/iA, oscillated due to quan- 
tum tunneling between clockwise and counterclockwise 
directions. The angular momentum associated with the 
current was of order L c ^ 10 10 7i. The total angular mo- 
mentum of the system is L to t = L c + L m , where L m is 
the mechanical angular momentum of the solid matrix 
bearing the current. To conserve L to t, oscillations of L c 
between L c — 10 10 ?i and L c = ~lO la h must be accompa- 



nied by simultaneous torsional oscillations of the solid 
matrix between the states with the angular momenta 
L m = -10 10 ?i and L m = I0 10 h, so that L tot = 0. One 
can argue that the current loop is always firmly attached 
to a substrate which is firmly attached to another solid, 
etc., so that the angular momentum is transferred to the 
infinite mass, like the linear momentum in Mossbauer 
experiment. If this were true, the SC current would 
simply bounce clastically between clockwise and coun- 
terclockwise directions, making conservation of the total 
angular momentum irrelevant. It is easy to see, how- 
ever, that for a current oscillating at a high frequency, 
e.g., fo ~ 10 9 — 10 10 s _1 , the analogy with the Mossbauer 
effect breaks down. Indeed, during the period of oscilla- 
tions, the elastic deformation cannot travel more than a 
distance v s / fo away from the current; v s being the speed 
of the transverse sound. For, e.g., w s ~5xl0 3 m/s and 
/o~5xl0 9 s~ 1 one obtains v s /fo ~ 1 /im. Thus, the part 
of the solid involved in the conservation of the angular 
momentum is small, not macroscopically large. Con- 
sequently, the torsional oscillations of the part of the 
solid matrix "co-wiggling" with the current must gener- 
ate phonons of frequency fo in the surrounding matter. 
As long as the speed of sound is finite this should result 
in a decoherence of quantum oscillations of the current. 
Similar effect exists due to the interaction of the magnetic 
moment of the current with vacuum photons. The differ- 
ence from the phonon problem is that at fo ~ 5x 10 9 s _1 
the wavelength of the light, A/ = c/ fo, is about five or- 
ders of magnitude larger than A s . Consequently, the 
vacuum properties of the photons depend strongly on the 
geometry of the experiment, in particular, on the metal 
shielding of the SQUID. 

The above picture is quasiclassical. From the quantum 
mechanical point of view, the states of the tunneling SC 
current are not classified by specific values of the an- 
gular momentum. The interaction of the current with 
the solid matrix entangles the angular momentum states 
of the current with the angular momentum states of the 
matrix. In the absence of the external noise, the total an- 
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gular momentum is, of course, a well-defined conserved 
quantum number. In this respect, it is important to un- 
derstand why the external noise acting on the macro- 
scopic solid matrix (like, e.g., vibration of the building) 
does not instantaneously destroy the quantum entangle- 
ment between the angular momentum states of the cur- 
rent and the angular momentum states of the matrix. We 
will discuss this issue in the context of the local versus 
global mechanisms of decoherence. 

The paper has the following structure. The decoher- 
ence due to the exchange of the angular momentum with 
the solid matrix is studied in Section II. The limits in 
which the dimensions of the SQUID loop are smaller 
and greater than the wavelength of the emitted phonons, 
X s = Vs/fo, are considered in the Subsections II-A and 
II-B correspondingly. The role of the external mechani- 
cal noise is discussed in the Subsection II-C. The photon 
effects are studied in Section III. The decoherence in the 
unshielded and shielded SQUID are computed in Subsec- 
tions III- A and III-B, respectively. Conclusions that may 
be useful for building SQUID-based qubits are summa- 
rized in Section IV. 



II. DECOHERENCE DUE TO PHONONS 

Consider quantum oscillations of the current J in a 
flat loop of spacial dimension R at a frequency / . The 
borderline between the cases of small and large SQUID 
studied in this Section is the dimension R ~ A s . For fo 
of the order of a few GHz, it is R ~ 1 /im. Obviously, 
both cases, R < A s and R > X s must be of interest 
for the ongoing experiments on quantum superposition of 
SQUID states and in connection with the goal of building 
a SC qubit. 

Let the magnitude of the dimensionless angular mo- 
mentum associated with the current be L. The magnetic 
moment of the current is 

M = [i B L = Ja/c , (1) 

where (1b is the Bohr magneton and a is the area of 
the loop. Thus, for a given current and given dimen- 
sions of the SQUID the value of the oscillating angu- 
lar momentum of the SC electrons can be found from 
L = Ja/fisc. For, e.g., a 0.1 /xm (small) SQUID carrying 
an electric current of J ~ 0.1 /iA, this gives L ~ 10 2 , while 
for a 100 /im (large) SQUID carrying an electric current 
of J — 10 fj,A, one obtains L ~ 10 10 . 

The dynamical torsional deformations of the lattice 
are described by the transversal displacement field u(r, t) 
satisfying 

V • u = . (2) 

These deformations do not affect the density of the ionic 
lattice and, in the long-wave limit, do not modify elec- 
tronic states. The minimal, allowed by symmetry, cou- 
pling of such deformations with the electric current j(r, t) 



must be proportional to J d 3 ryu. Let us show that the 
coefficient of proportionality is m e /e, where m e and e are 
bare electron mass and charge. Without loss of general- 
ity, one can employ the following argument in order to 
derive the interaction term. Consider the classical den- 
sity of the kinetic energy of electrons and ions, 

KE = n e — — ^ + m — - — , (3) 

where Tig A 3X6 concentrations of electrons and ions, Mj 
is the ionic mass, and v e is the velocity field of electrons 
in the laboratory coordinate frame. Notice now that the 
electronic current is formed by the electron band states 
in the coordinate frame co-moving with the lattice. Con- 
sequently, 

j = en e (v e - u) . (4) 

Expressing Eq. (|^) in terms of physical variables j and u, 
one ends up with the coupling of the form (m e /e)(j-u). 
The generalization of this argument for quantum opera- 
tors is trivial. It gives the following Hamiltonian of the 
effective interaction in the laboratory frame: 

m f 

n int = <?ryu ■ (5) 

This formula is a consequence of the translational sym- 
metry and is, therefore, parameter-free. It can be used 
as long as Timt results in a small perturbation of one- 
electron energy band. 

A. Small SQUID 

We shall start with the case of a small current loop 
of R < A s , in the XY plane. The shape of the loop 
will be irrelevant. In the context of its magnetoelastic 
interaction with the solid matrix it is equivalent to a 
magnetic atom of angular momentum L, imbedded in a 
solid. We shall denote the states with L up and L down 
by | T ) an d | | ) respectively. They are the eigenvalues of 

L,|T) = i|T>, L z \i) = -L\i). (6) 

With no regard for the conservation of the angular mo- 
mentum, the ground state of the system can be approxi- 
mated by 

|o) = -±=(|T) + |4», (7) 

while the first excited state, with the same accuracy, is 

|i) = ^(IT)-U)). (8) 

Let the energy separation of these two states be A — 
hujQ = hfo. If the energy scale of the experiment does 
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not significantly exceed A, then any quantum state of 
the SQUID is a superposition = Ci|0) +C 2 |1). One 
of the mechanisms of decoherence of this superposition, 
that exists down to zero temperature, is the decay of |1 ) 
onto |0) accompanied by the radiation of the quantum 
of energy A. 

As has been discussed in the Introduction, the conser- 
vation law requires that the oscillations of the current 
are accompanied by the torsional oscillations of the solid 
matrix, so that the total angular momentum stays con- 
stant (e.g., zero). Such a local wiggling of the matrix 
must result in the finite probability of the emission of a 
transverse phonon of frequency /o- Since the wavelength 
of the phonon, X s , is large compared to the dimensions of 
the SQUID, its effect on the SQUID is equivalent to the 
uniform local .rotation of the solid matrix at the position 
of the SQUIDEJ. In terms of the deformation field u(r, i), 
the angular velocity of this rotation is given by 



i. 



-Vxu 



(9) 



It is important to notice that the decoherence rate of 
Eq. (0) is proportional to the fifth power of / and to 
the fourth power of the size of the SQUID. It is inversely 
proportional to the fifth power of the speed of the trans- 
verse sound, making it important to use solid matrices 
(substrates) of high shear modulus. For practical values 
of the parameters: k B T<A, R^O.lpm, J^O.lpA, 
p^5g/cm 3 , v s ~5xl0 3 m/s and /o~5xl0 9 s _1 (that 
is, A s ~l/im), Eq. (14) gives r s ~ls _1 . Consequently, 
the above mechanism of decoherence should not be of 
great concern for a small SQUID. For a large SQUID the 
situation will be quite different, as is discussed below. 

We should emphasize that Eq. ( |T4| ) can be used for the 
estimate of the decoherence rate only at R < X s =v s / Jo- 
in experiments of Ref. || and Ref. || the size the SQUID 
was large in comparison with A s . The interaction of 
such a SQUID with phonons is non-local and cannot be 
treated by the above method. The non-local theory of 
the phonon emission by a large SQUID is developed in 
the next Subsection. 



Accordingly, the lattice velocity field at the position of 
the SQUID is u = f2 x r. Substituting this into Eq.(||), 
one finds 



n 



where 



eff 



r r x j 



(10) 



(11) 



stands for the angular momentum of the SC current. The 
effective interaction (|l^), is mandated by symmetry and 
is, therefore, parameter free. Correspondingly, the mech- 
anism of decoherence provided by this effect is universal. 

Based upon Eq. ( |l0| ) the rate of the transition from |1 ) 
to |0 ) is given by 



r s = 2h < 0|L|1 > [J e ™(A)l < 1|L|0 > 



(12) 



where J env (A) is the spectral function of the environ- 
mental coupling for phonons, 

J env {A) =7r5^<k I *|n|0)(0|n|k J t)*(A - to.*,-) ; (13) 



k, i, and wti denote the wave vector, polarization, and 
the frequency of the phonon, respectively. Further com- 
putation along the lines of Ref. |ll| yields 



KT 2 

r s (T) = -_fc5 C oth 



2k B T 



(14) 



J 2 a 2 



Sirfre 2 



K coth 



A 



2k B T 



where p is the mass density of the solid matrix and k s — 
2tt / X s is the wave number of the emitted sound. 



B. Large SQUID 

For R > X s , one should employ the general form of 
the interaction given by Eq. (|^). This term results in the 
coupled dynamics of the currents and the lattice displace- 
ments. However, in reality, the large difference between 
the ionic and electron masses makes the renormalization 
the SC dynamics insignificant. In what follows, we will 
ignore the effect of Eq. (^) on the spatial and temporal 
structure of the current and the SQUID flux <E> gener- 
ated by the current. We shall be concerned with the fact 
that, due to Eq. (||), the currents serve as the source of 
phonons in the elastic equation. In other words, the solid 
lattice must take the recoil from the oscillating current. 
This effect is mandated by conservation laws and it leads 
to the decoherence of the quantum dynamics of the flux. 

Inside a good conductor, either metal or a supercon- 
ductor, the longitudinal electric fields are screened with 
a typical time scale of the plasma oscillations, t ~ 10~ 15 s. 
The flux dynamics is much slower. Consequently, the lon- 
gitudinal phonons that change the local concentration of 
ions should be excluded from our consideration. This can 
be done by supplementing Eq. (@) with the condition (j|) . 
For the purpose of estimates we shall adopt the simplest 
model of uniform and isotropic elastic medium. Then, 
the energy of the free transverse phonon field is 



T~Cph 



d 3 ? 



-pu 



(15) 



where p is the shear modulus of the solid and U{j = 
^(diUj + djUi) is the strain tensor (^ i uu being zero for 
transverse phonons). In this model the torsional strains 
are described by just one elastic modulus p = pv 2 . Ac- 
cordingly, the transverse sound velocity v s is indepen- 
dent of the phonon polarization (which is orthogonal to 
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the phonon wave vector). We shall further simplify our 
consideration by neglecting all differences in the actual 
material composition of the experimental setup, that is, 
by assuming that the phonon spectrum is the same inside 
and outside the part of the solid matrix that carries the 
SC current. This assumption, while not valid in experi- 
ment, should not significantly affect our estimate of the 
decoherence rate. 

The canonical quantization of the phonon field yields: 




(16) 

Here II = pii is the momentum of the phonons that is 
canonically conjugate to u, uj^i — v s k is the frequency of 
the phonon of the wave vector k and polarization i, and 
V is the volume of the system. Due to the isotropy, Wki 
for the transverse phonons (k-e, = 0) does not depend 
on the polarization. Substituting Eqs. into Eq. ( |l5| ) 
and Eq. (§) one obtains 



n = n 



ph 



+ 



E 



/ fruj 



1/2 



(jk-e 4 )(a kl -ak) 



(17) 



where jk = / d 3 r j exp(— ikr) is the spatial Fourier com- 
ponent of the current density j (r) , and the summation is 
over k and i satisfying ke^ = 0. The Fermi golden rule, 
then, yields the following expression for the decoherence 
rate: 

2tt m 2 
T ph = — t-^- A coth 



h 2e 2 p 



E 



d 3 k 



2k B T 



KOKjk-eOli)! 2 ^-^) 



(18) 



where | 0) and 1 1 ) are given by Eqs. (0) and (||). 

We shall now compute jk- The SC current can be writ- 
ten in terms of the SC phase (p and the vector potential 
A, 



47TA! 



Z7T 



(19) 



with $o = hc/2e and Xl being the flux quantum and 
the London penetration length, respectively. The vector 
potential satisfies the Maxwell equation, 



«-, «-, » 47T. 

V x V x A = — j , 

c 



(20) 



with the current density j given by Eq. (|T^) inside the 
SC loop and by j = outside the loop. 

To simplify calculations we shall adopt the ring geom- 
etry of the SQUID. In the presence of the flux the 
phase ip winds around the ring by (pj = 27r$/$Oj which 
is the Josephson phase in the junction cutting the ring. 
We shall study the problem in cylindrical coordinates 
[z, r, 4>), with the Z-axis passing through the center of 
the ring perpendicular to its plane, r standing for the ra- 
dial coordinate, and <fi being the polar angle in the plane 
of the ring. The solution for the phase is 



VJ<f> 
2tt 



$ ' 



(21) 



Due to the cylindrical symmetry of the currents, the 
only non-zero component of the vector potential is A = 
A^r, z). Then, inside the ring, Eq. ( pp| ) reduces to 



d_ 

dr 



1 d , „x 

-jr(rA) 
r or 



d 2 A 

d 2 z 



+ X7 2 A 



$ 1 



2TT\j 



while outside the ring one has 



d_ 

dr 



1 9 , AS 

r or 



2 A 

d 2 z 







(22) 



(23) 



These equations must be accompanied by the boundary 
conditions for A and for non-zero components of the mag- 
netic field, 



H r = - 



dA 

dz 



H z = 



1 d(rA) 
r dr 



(24) 



In this paper we shall not pursue the exact solution of the 
problem for the finite cross-section of the ring carrying 
the current. Instead, we will make use of a thin-ring ap- 
proximation in which the thickness of the ring yb, where 
b stands for the area of the wire crossection, is small com- 
pared to its radius R as well as to Xl ■ Then, in cylindrical 
coordinates, the only non-zero component of the current 
density is j<t,(z, r). It equals J/b inside the ring and zero 
outside the ring. At k\fb -C 1 the Fourier transform of 
such a distribution of the current is 



jk = —i2nRJi(kj^R)Jn k 



(25) 



where rik J- k is the unit vector in the plane of the ring, 
k± = fcsinff, 8 is the angle between k and the Z-axis, 
and Ji(k±R) stands for the Bessel function of the first 
order. 

The quantization procedure consists of assigning the 
operator J to the total current J. In the two-level ap- 
proximation, one introduces the states |±) of the cur- 
rent operator such that ±J are the respective eigenval- 
ues: J|±) = ±J|±). In terms of the angular momentum 
operator, these states are identical to those in Eq. (^), 
that is |+) = | t) an d |— ) = I I)- Tunneling between 
these two degenerate states produces new states, Eq. (Q) 
and Eq. (0), which are split by the energy A. These 
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states are characterized by the zero current, ( |J| ) =0 
and (1|J|1) = 0. The transition matrix element is 
(0|J|1) = J. Substituting Eq. @ into Eq. @ we 
get 



2™ 2 J 2 R 2 
1 = — ^ coth 



Tie 2 



fcf / dcos6 J 2 (k s Rsint 
Jo 



(26) 



The limit of /c s i? <C 1 corresponds to a small SQUID. It 
is easy to see that in this limit T p h is proportional to kg 
and Eq. ( p6| ) becomes Eq. (|l4|). In the limit of a large 
SQUID, k s R > 1, we find 



ph 



Trm 2 J 2 R 
he 2 p 



k 2 coth 



A 



2k B T 



(27) 



Based upon Eq. (]27|), let us make an estimate of the 
decoherence rate for, e.g., the experiment of Ref. H. 
At k B T<A, for R~ 0.1 mm, J~3/xA, p~8g/cm", 
w s ~ 5xl0 3 m/s and / ~ 2xl0 9 s^ 1 (that is, A s ~ 2.5 /urn), 
Eq. ( p7j ) gives I 1 },/, ~ 10 6 s -1 . This is a significant deco- 
herence rate that would limit the quality factor of the 
corresponding qubit by the value of about one thousand. 



C. Global noise 

Here we will compare the effect of the global noise — 
uncontrolled rotations of the solid matrix as a whole at 
some angular velocity ^g(^) — with the above estimates 
for the local effects due to phonons. The Z-component 
of the global rotation removes the degeneracy between 
clockwise and counterclockwise current states. This is a 
particular case of the Barnett effect: A rotating solid de- 
velops magnetisation proportional to the angular velocity 
of the rotatiorJl3. In application to the SQUID this effect 
can be described by a two-state Hamiltonian written in 
the rotating coordinate frame: 



H G = -&s x -2hm G (t)s z 



(28) 



where s XjZ are spin-1/2 operators. We want to estimate 



the effect of the second term in Eq. (E8|) on coherent os- 
cillations of the SC current. (Notice that for externally 
imposed rotations, the sign of this term is opposite to the 
sign of Eq. (|^|) that was written for phonons dynami- 
cally produced by SQUID oscillations in the laboratory 
frame.) For a macroscopic solid matrix, the characteris- 
tic correlation time of (f2G(i)f2£?(0)) cann °t be less than 
the time it takes the sound to travel across the matrix. 
Consequently, on the time scale £~l//o, random rota- 
tions of the equipment as a whole are slow enough to 
permit the treatment of Eq. ( p8| ) within the adiabatic 
approximation. 

In the adiabatic approximation the eigenvalues and 
eigenf unctions of the Hamiltonian (E8I) are 



e± = ±VA 2 + [2hm G ] 2 
^± = ^[c q =(t)|t)±c±u; 



V2 

where C± are given by 



'1 ± 



2hm 



G 



y/A 2 + [2hm c 



(29) 



(30) 



One can see that random rotations of the system as a 
whole do not significantly perturb the states | ) and | 1 ) 
given by Eq. (0) and Eq. (|) only if VL G satisfies 



(31) 



For a large SQUID with L ~ 10 10 and uj ~ 10 10 s _1 , this 
gives rig ^ 1 s ; which must be of practical importance 
in the situation when the equipment is subjected to ran- 
dom movements. For a small SQUID, with small L and 
practical values of cjo, the condition ( |3l| ) can be satisfied 
by a very large margin. 

The time dependence of il G is another factor. The 
adiabaticity implies that 



A 



4LIO 



> 1 



(32) 



G 



Then, the main effect of the time dependence of the 
global rotations is generation of additional harmonics in 
the Rabi oscillations between clockwise and counterclock- 
wise SC currents due to the time dependence of e± . The 
corresponding decoherence rate, Tq, can be estimated 
from the variation of e± during one cycle, 2tt/ujq, of the 
undisturbed Rabi oscillations. The global random rota- 
tions occur when the equipment is subjected to a random 
external torque. Let this torque result in an angular ac- 
celeration a = Sic, so that the change of £Iq during one 
cycle is Sfla ~ a/tug. This gives 



T G ~5e + /h~u a /A 



(33) 



where A — tilo 2 ,/ AL\a\. Thus, the quality factor for the 
above mechanism is Q = A. It is entirely determined by 
the adiabaticity of global rotations. 

If the typical frequency of the mechanical noise is ui n , 
then a^uj n D,G- Assuming that the condition ( f3l| ) is sat- 
isfied, this implies 



Q = A > LU / CJ„ 



(34) 



In principle, at small ujq, the effect of the mechanical 
noise can overpower the decoherence from local phonon 
effects, which decrease as some power of u>q. However, 
for any practical values of luq and a in the coherence 
experiments, the adiabadicity factor A is very large, so 
that the decohering effect of the external global noise can 
be safely ignored. 

We also would like to make the following interesting 
observation. Consider rotations flc = at that last long 



5 



enough to violate the condition ([H]). According to Eq. 
(^p|), as time goes from — oo to +00, the states tp± of 
Eq. ( p^ ) switch between | | ) and | J,). This, however, is 
true only in the limit of \a\ — >0. At finite \a\ the answer 
depends on the adiabaticity factor A. If at t = —00 the 
SQUID is prepared in, e.g., the state with the clockwise 
SC current, then, the probability for the SQUID to switch 
to the counterclockwise cuixent at t = +00 is given by 
the Landau- Zener formulalij 



P, 



LZ 



1 



(35) 



This probability is high if the angular acceleration sat- 
isfies a < a c = ttujq/AL. Consequently, at low tunnel- 
ing rate the relatively slow mechanical rotation can pro- 
vide the quantum-mechanical switching between clock- 
wise and counterclockwise currents. It should be noted, 
however, that the Landau-Zener transitions can also be 
generated by the magnetic field. Thus, the above ef- 
fects of uniform rotations can only be observed if the 
SQUID is shielded from the magnetic fields with an ac- 
curacy hbH < %£Iq. 



III. DECOHERENCE DUE TO PHOTONS 

The problem of decoherence due to the emission of 
photons of frequency /o is very similar to the problem of 
the emission of phonons. The main difference is that the 
vacuum wavelength of the light, A; = c//o, is typically 
large compared to the size of the SQUID that exhibits 
quantum oscillations of the current. The electromagnetic 
radiation by such a SQUID into the open space is equiv- 
alent to the radiation of a point magnetic dipole. If, 
however, the SQUID is shielded by a metal placed at a 
distance that is comparable to or smaller than A;, the 
decoherence rate becomes strongly geometry-dependent. 
These two problems are considered in the following two 
Subsections. 



A. Decoherence in the open space 

If the wavelength of the emitted photons is large com- 
pared with the SQUID size, then the SQUID can be 
treated as a point particle with an angular momentum L 
which is perpendicular to the SQUID loop. This angular 
momentum interacts with the photon field via Zceman 
Hamiltonian 



n 2 



-^bL z H z 



(36) 



where H z is the Z-component of the magnetic field of 
the vacuum photons. We are interested in the transition 
between the tunnel-splitted quantum states given by Eq. 
(|) and Eq. (§). The expression for the rate is similar to 
Eq. (|T2|), where the spectral density of the photons can 
be obtained by either quantizing H(r) or taking 
from the theory of electromagnetic fluctuations! 14 ''' 1 



The Fermi golden rule then yields the following expres- 
sion for the decoherence rate 



r*(T) 



— n%L z kf coth 
^LjVfcfcoth 



2k B T 
A 



(37) 



2k B T 



where fc; = A/ he is the wave vector of the emitted light. 
It is proportional to the third power of /o as compared 
to the fifth power of fo in Eq. (|l4]). The reason for the 
difference is the additional time derivative of the boson 
field in the Hamiltonian of Eq. (hCl) as compared to Eq. 



3- 

Notice that the identical result for Y z follows from 



r*(o) 



i 

Hloq 



wherel 







(38) 



(39) 



is the intensity of the magnetic dipole radiation due to the 
classical dynamics of the magnetic moment M = /i^L. 

For a small SQUID (as defined above in Sec. II) 
with i2~A s ~100nm and J-O.l^A, Eq. ((371), at 



/o~10 s and k B T < A, gives negligible decoher- 
ence, r.z ^ 1CU 8 s -1 . However, for a large SQUID with 
fl~100/im and J~3^A, at /o^lO^s" 1 and k B T < 
A, one obtains rz~10 7 s -1 . Thus, for a large SQUID, 
the radiation of photons into the open space can easily 
reduce the quality factor of the SQUID down to one hun- 
dred. 

The decoherence rate can be decreased by choosing a 
double- loop geometry with equal areas of the single loops 
and equal currents flowing in the opposite directions, as 
was actually done in Ref. a In that case, the total mag- 
netic moment of the system is zero and the radiation is of 
the quadrupolar nature. If the magnetic moment is com- 
pensated exactly (which must be difficult to achieve in 
experiment) the radiation rate will be reduced by a factor 
(kR) 2 < 1. At a finite compensation, 7 = AM/M < 1, 
the decoherence rate ( |37| ) acquires a factor r ) 2 . 

Here we have neglected the effects of dc and low- 
frequency adiabatic ac magnetic fields H{t). These ef- 
fects are equivalent to the effects of global rotations stud- 
ied in the Subsection II-C. To estimate them quantita- 
tively one should replace Qq and a by lo b = eH/2m e c 
and Lo ac LO B respectively, where oj ac is the typical fre- 
quency of the ac field. Then, conditions ( J3l| ) and ( p^ ) 
become 



and 



A 



ALuJacUH 



» 1 



(40) 



(41) 



G 



respectively. In a two-loop design, where the magnetic 
moment is compensated by a factor 7 = AM/M < 1, 
L in Eqs. (|o|) and (|l]) should be replaced by jL. 

For, e.g, a small SQUID with L ~ 100 and no com- 
pensation, Eq. (El) at luq ~ 10 10 s _1 translates into <C 
10 Oe. For a large SQUID with 7 = 10~ 2 , L = 10 10 and 
ujq ~ 10 s _1 , the fields that do not disturb the states | ) 
and I 1 ) should satisfy H < 10~ 5 Oe. If Eq. © is sat- 
isfied, then substituting it into Eq. ( f4l| ) one obtains the 
following relation for the quality factor coming from the 
low-frequency magnetic noise alone: 



Q = A > UJo/iOa 



(42) 



B. Decoherence in the presence of metal shielding 

We shall now study the case when the SQUID loop is 
adjacent to a metal sheet parallel to the plane of the loop. 
Following published experiments, we shall assume that 
the distance between the SQUID and the shielding, d, is 
much smaller than the wavelength of the vacuum elec- 
tromagnetic radiation A/. In that case the radiation be- 
comes strongly renormalized by the conducting medium 
and the formulas of the previous Subsection can no longer 
be used. Now the main source of decoherence is the dis- 
sipative current in the metal shielding induced by the ac 
fields of the SQUID. Correspondingly, the decoherence 
rate can be computed as 



M 



P 



(43) 



where P is the power absorbed by the shielding. Let 
H = H w (r) exp(iwo^) be the magnetic field generated by 
the_oscillating current in the SQUID. Then P is given 

b y y 



P 



16tt V 2tr7 



|H w | 2 )d/, 



(44) 



where a is the electric conductivity of the shielding, (...) 
means quantum-mechanical average, and the integration 
goes over the metal surface facing the SQUID. 

Eq. ( phi] ) can be used when the thickness of the skin 
layer, 8 — c/ \/2tt(ju)q, is small compared to the thick- 
ness of the shielding metal, D, but large compared to 
the mean free path of electrons of the metal, Iq. These 
conditions were apparently fulfilled in the experiment of 
Rcf. H for which we estimate 6 ~ 1 /im and Iq < 0.1 /im 
at D ~ 8 /im. The condition d <C .Aj allows one to use 
the quasistationary approximationEj to obtain H. In 
this approximation the field is formed by two current 
loops, one being the mirror image of the other with re- 
spect to the surface of the shielding. At the metal sur- 
face, z = d, this field has the tangential component only, 
H r (r) . For a thin circular loop carrying the electric cur- 
rent J(t) = J exp(iwoi), it is given by 



H r = 



4JRd 



COS ( 



o (d 2 + r 2 + R 2 - 2rR cos (j>) 3 / 2 ' 



(45) 



which can be expressed in terms of the elliptic integrals^ 
Substituting Eq. (E|) into Eq. (EJ), one obtains 



2 J 2 



hc V°7o 

where the function F is given by 

/•oo r p-k 



F(R/d) 



(46) 



+ x 2 + £ 2 - 2x£ cos 0) 3 / 2 



(47) 



In the two limiting cases of small and large d, the geo- 
metrical factor (WTu reduces to 



F = — , d<i? 
4a 



(48) 



F 



3a 2 
32d4 



d> R 



Here we have introduced the length and the area of the 
current loop, / = 2irR and a = irR 2 , respectively, in 
order to emphasize the fact that the above limiting ex- 
pressions are correct for flat current loops of arbitrary 
shape. Numerical analysis shows that these expressions 
hold for d/R < 0.3 and d/R > 5, correspondingly. For 
the double-loop with a compensated magnetic moment 
due to equal single-loop currents flowing in the oppo- 
site directions, the decoherence rate practically does not 
change in the limit of d <C R. In the opposite limit of 
d 3> R, it reduces by a factor (R/d) 2 . 

Note that the frequency dependence of the rate ( |4^ ) 
follows from the frequency dependence of the skin depth, 
<5 l/\/7o- Then the Maxwell equation gives E cx 
8-(dH/dt) cx y/fo for the electric field in the skin layer 
and the dissipation rate due to Joule's power losses, 
<jE 2 5, becomes proportional to \fja- Divided by /q in 
Eq. (f43j), it gives the l/VTo dependence of the decoher- 
ence rate. If the thickness of the metal shielding D is 
smaller than 6, the electric and magnetic fields are not 
significantly modified by the shielding, so that E cx foH , 
with the proportionality factor determined by the geome- 
try of the SQUID and its distance to the shielding. Thus, 
in the low- frequency limit (D < 6), the dissipation in the 
shielding P~aE 2 D is proportional to f 2 and Tm due to 
shielding is proportional to /q. 



/o~2xlO a 



and 



For J — 3 uA, a 
R^d, Eq. ( flo] ) gives Tm ~ 10 9 s _1 , which was, probably, 
the case in the experiment of Ref. ||. This shows that for a 
SQUID carrying a microampere current the above mech- 
anism can provide a very high decoherence rate. Notice 
that Tm can be drastically reduced by increasing the dis- 
tance, d, between the SQUID and the shielding. Indeed, 
according to Eqs. @, Tm <x d~* at R < d < A;. Too 
large d, however, would reduce the effectiveness of the 
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shielding in protecting the SQUID from external radio 
signals. Choosing smaller SQUIDS operating at smaller 
currents should be more beneficial for qubit designs. 



IV. CONCLUSIONS 

We have studied generic mechanisms of decoherence 
mandated by the conservation laws - - emission of 
phonons and photons of the oscillation frequency fo- Our 
practical conclusions are as follows: 

The decoherence due to the above mechanisms scales 
as the second power of the current. 

For small SQUIDS of size R < v a / fo, the decoherence 
due to the emission of phonons at T < A is negligible. 

For large SQUIDs of size R » v s /f , the emission of 
phonons can significantly limit the quality factor. The 
corresponding decoherence rate scales linearly with the 
size of the SQUID and quadratically on the oscillation 
frequency. 

In the absence of the metal shielding, the emission of 
photons is negligible for small SQUIDs but becomes sig- 
nificant for large SQUIDs. It scales as the fourth power 
of the size of the SQUID and as the third power of the 
oscillation frequency. 

Decoherence due to the shielding strongly depends on 
the geometry of the experimental setup. It may com- 
pletely destroy the coherence in large SQUIDs and can 
be the main mechanism of decoherence in small SQUIDs. 
The shielding must be provided by a metal sheet of thick- 
ness greater than the skin layer at the oscillation fre- 
quency fo- To achieve small decoherence, the distance 
to the shielding, while small in comparison with c/ fo, 
should be considerably greater than the loop size. 

The effects of the external mechanical and magnetic 
noise are proportional to the total magnetic moment of 
the SQUID, making small SQUIDS less susceptible to the 
noise than large SQUIDs. 
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